Functionnalized probes for AFM force spectroscopy: eigenmodes shape and stiffness 
calibration through thermal noise measurements 
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The functionalization of an Atomic Force Microscope (AFM) cantilever with a colloidal bead 
is a widely used technique as soon as the geometry between the probe and the sample has to be 
controlled, in particular in force spectroscopy. But some questions remain: how does a bead glued 
at the end of a cantilever influence its mechanical answer? And more important for quantitative 
measurements, can we still determine the stiffness of the AFM probe with traditional techniques? 

In this article, the influence of a colloidal mass loading on the eigenmodes shape and resonant 
frequency is investigated by measuring the thermal noise on rectangular AFM microcantilevers with 
and without a bead attached at their extremity. The experiments are performed with a home-made 
ultra-sensitive AFM, based on differential interferometry. The focused beam from the interferometer 
probes the cantilever at different positions and the spatial shapes of the modes are determined 
up to the fifth resonance, without external excitation. The results clearly demonstrate that the 
first eigenmode almost doesn't change by mass loading. However the oscillation behavior of higher 
resonances present a marked difference: with a particle glued at its extremity, the nodes of the mode 
are displaced towards the free end of the cantilever. These results are compared to an analytical 
model taking into account the mass and the inertial moment of the load in an Euler-Bernoulli 
framework, where the normalization of the eigenmodes is explicitly worked out in order to allow a 
quantitative prediction of the thermal noise amplitude of each mode. A good agreement between 
the experimental results and the analytical model is demonstrated, allowing a clean calibration of 
the probe stiffness. 
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I. INTRODUCTION 



Atomic Force Microscopy (AFM) is currently used in 
a great variety of studies from various disciplines to 
measure small forces by measuring the deflection of a 
microcantilevcr [1]. In biophysics for example, it has 
been applied to the unfolding of protein [2, 3], prob- 
ing the structure of biological membranes [4], monitor- 
ing the mechanical response of living cells [5, 6]... In 
nanotechnology as well, micro-scale levers find applica- 
tions in Micro-Electro-Mechanical Systems (MEMS) and 
other nanotechnological devices [7]. In material, surface 
or nano sciences in general, AFM probes appear as a 
cornerstone for quantitative studies at nanoscale [8] . 

All those applications exploit the great accuracy in 
measuring the cantilever deflection offered by AFM and 
converting this measurement in units of force assuming 
the cantilever behaves like a spring with known stiffness. 
Manufacturers often specify the spring constant of their 
cantilevers in a wide range of values, mainly because of 
the great uncertainties in the dimensions, particularly 
the thickness, resulting from the fabrication process. To 
overcome this problem several techniques have been pro- 
posed to calibrate the cantilever spring constant [1, 9-12]. 
The reader is referred to the work of Burnham and co- 
workers [13] and the references therein for a comparative 
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summary of the different techniques. 

One of the first and still most commonly used calibra- 
tion method is the so-called thermal calibration method 
based on the measurement of the vibration amplitude of 
the free end of a cantilever exited by thermal noise [11]. 
The first peak of the thermal noise spectrum is related 
back to the spring constant of the cantilever modeled 
as an harmonic oscillator. In a more accurate model, 
Butt and Jaschke [14] introduced a correction factor de- 
duced from the Euler-Bernoulli description of the flexu- 
ral dynamic of a free-clamped beam. In a previous work, 
we demonstrated that measuring this thermal noise for 
the first resonant modes of the cantilever provides an 
excellent benchmark to probe the mechanical response 
of the cantilever and compare it to a simple mechanical 
model [15]. 

In this article, we extend this method to the case of 
functionalized AFM cantilevers. Indeed, as soon as the 
geometry between the probe and the sample must be con- 
trolled, it is common to use a colloidal bead fixed at the 
free end of the lever. The radius of curvature of the "tip" 
is then controlled and stable, and offers a clean sphere- 
plane geometry to study interaction at nanometric dis- 
tances. These modified probes (whether home made or 
now commercially available) arc commonly used in force 
spectroscopy, in particular for the measurement of the 
nanorheology of confined fluids [16, 17] or of the Casimir 
interaction [18]. How does this loaded mass influence the 
mechanical response of the AFM cantilever? Can the 
common techniques (thermal noise calibration in partic- 
ular) to determine the stiffness still be used? 
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In this work, we measure the thermal noise spectra of 
the cantilever deflection on its whole surface and com- 
pare the rms amplitudes obtained with and without a 
bead loaded at its free extremity. Furthermore, we com- 
pare the results to a simple mass model which modifies 
only the boundary conditions of the classical rectangular 
beam theory [19]. A good agreement between experimen- 
tal data and this analytical model will be demonstrated, 
showing that thermal noise calibration of the probe stiff- 
ness is still perfectly pertinent for such cantilevers. 

The paper is organized as follows. Section II describes 
the theoretical approach, with a special emphasis on 
eigenmodes normalization to allow the prediction of ther- 
mal noise amplitude of each mode. Section III details the 
experimental results in the light of this model, for two 
cantilevers probing various mass ratios between the can- 
tilever and the colloidal bead. Eventually, a discussion 
and conclusion is given in section IV, with a specific fo- 
cus on how our conclusions can be applied to the classic 
angular deflection measurement technique. 

II. ANALYTICAL DESCRIPTION OF 
THERMAL NOISE 
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A. Flexural eigenmodes of a clamped cantilever 

In the Euler-Bcrnoulli framework to describe the mi- 
crometer sized mechanical beam, we assume that the can- 
tilever length L is much larger than its width W, which it- 
self is much larger than its thickness T (see figure 1). The 
flexural modes of the cantilever are supposed to be only 
perpendicular to its length and uniform across its width. 
The deformations can thus be described by the deflec- 
tion Z(x, t), with t the time, and x the spatial coordinate 
along the beam normalized to its length L. The bead is 
supposed to be non deformable, and thus described as a 
rigid mass mj, at the free end of the cantilever. It may 
be offset from the neutral axis, and carrying a non negli- 
gible inertia: these effects are taken into account with an 
inertial moment m&? "jf, r g being the equivalent gyration 
radius computed at the free end of the ca ntile ver. For a 
sphere or radius r, r g is then equal to yf/5r. We ne- 
glect in our analysis any offset along the cantilever axis, 
and any coupling with torsion. Figure 1 sketches the ex- 
periment and the applied model. Following [19], we will 
include the effect of the bead in the boundary conditions 
of the cantilever dynamics. 

The equation of motion for the cantilever, once the 
variables in time and space separated, can be written 



k d 4 z 
3d? 



m c uj z 



(1) 



with k the static stiffness of the cantilever, m c its mass, 
Z(x,t) = z(x)e lu)t the deflection, and uj the pulsation. 
This equation can be rewritten as 



FIG. 1: Schematics of the experiment and analytical model. 
The cantilever (length L, width W , thickness H) is modeled in 
the Euler-Bernoulli framework by its neutral axis subject to a 
deflection Z(x, t). The bead glued at its free end is modeled as 
a point mass mj, located at the extremity. r g is the equivalent 
gyration radius of bead and accounts for the inertia in rotation 
of the bead at the cantilever end. In the experiment, the 
deflection can be measured at any position x and y (along 
and transverse to the cantilever axis) with a very low noise 
differential interferometer [15, 21, 22], sensing the optical path 
difference between the two laser beams represented on the 
figure. 



where is the spatial derivative of order n, and a is 
given by the dispersion relation: 



3m c w 2 
k 



(3) 



..(4) 



4 

a z 



(2) 



The generic solution to this equation is 

z(x) — acos(ax)+bsm(ax)— ccosh(aa;)— dsmh(ax) (4) 

The boundary conditions in x = corresponds to a 
clamped end, implying z(0) = ct zW(0) = 0, hence 
a = c and b = d. Defining R = b/a, the generic solution 
4 can thus be written 

z(x) = a (cos(acc) — cosh(o!ar) + R [sin(aa;) — sinh(ax)]) 

(5) 

In x = 1, the conditions on the force and torque are 
linked to the inertia in translation and rotation of the 
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bead [19]: 



z (3) (l) = -crW(l) 
z( 2 '(l)=a 4 mf 2 zW(l) 



(6) 
(7) 



where m = jm c is the mass of the bead normalized 



to that of the cantilever 



and : 



,/ L is the gyration 



radius of the bead r g normalized to the cantilever length 
L. Expressing those boundary conditions with expression 
5 leads to 



R = 



sin a — sinh a + am(cos a — cosh a) 
cos a + cosh a — am (sin a — sinh a) 
cos a + cosh a — a 3 fhf 2 (sin a + sinh a) 
sin a + sinh a + a 3 fhf 2 (cos a — cosh a) 



(8) 
(9) 



The values of a allowing this equality are quantified, 
and correspond to the spatial eigenvalues a„(?fi, f) of the 
resonant modes of the cantilever. They can be numeri- 
cally computed. For fh = (no bead), the last equation 
simplifies to the usual condition 1+cos a cosh a = 0, lead- 
ing to the common tabulated eigenvalues of a clampcd- 
free Euler-Bcrnoulli mechanical beam. The a n (rh, f) val- 
ues are reported for the fisrt 5 modes, for < fh < 2 
and < f < 0.1, in tables IV to VIII in appendix A. 
The corresponding shapes of the eigenmodes are plotted 
in figure 2. 

The length of the cantilever L, used in the normaliza- 
tion of x and thus impacting the spatial eigenvalues a n , 
is sometimes experimentally ill-defined due to the trian- 
gular shape of the cantilever end. The direct comparison 
of the experimental values of a n with the theoretical ones 
is thus hampered by this incertitude. However, their ra- 
tio is exempted from this bias, and can be used to check 
analytical predictions. In appendix A, we plot such ra- 
tios, useful to extract the values of fh and r from the 
experimental observations. 



B. Orthogonality relations and normalization 

The eigenmodes z n (x) are given by equation 5, with a 
dependance of the three parameters a, R and a on mode 
number n and on fh and f (the dependance in fh and f 
will be implicit in our notations): 

z n (x) = a n ( n (x) (10) 
= a n (cos(a„a;) — cosh(a n x) 

+ Rn [sin(a„x) — sinh(a„x)] ) (11) 

However, if we only consider this expression, then the 
orthogonality between two modes z n and z m does not 
hold: it is easy to show that for n =/= m, 

/ z n (x)z m (x)dx = -mz n (l)z m (l) - fhf 2 z^ (Vjz^ {I) 
Jo 

(12) 
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FIG. 2: We plot in this figure the 5 first normal modes (j>n{x) 
for various fh and f. The amplitude of the first mode decreases 
and the nodes of the higher order modes get closer to the free 
end when fh or f increase. 



Let us first focus on the case f 
the normal modes by 



0. We can then define 



(f> n (x) = z n (x) + z n (l)y / m5(x - 1) 



(13) 



where 8{x— 1) is the Dirac distribution centered in x = 1. 
With such definition, it is straightforward to prove the or- 
thogonality of the 4> n basis. Moreover, as 4> n (x) = z n (x) 
for any x € [0,1[, 4> n {x) thus obeys the initial diffcren- 
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tial equation 2 and the boundary conditions in the limit 
x = and x — > 1. Eventually, the normalization of </>„ is 
easily obtained by imposing 



<t>l(x)dx = a\ 



300 



C(x)dx + mC(l) =1 (14) 



where Cn(^) is defined in equation 10. This last equation 
thus imposes the values of a n to construct the orthonor- 
mal basis <j) n . The result of this process is illustrated in 
figure 2 for a few values of rh. 

Let us now give an energetic meaning to the Dirac 
term. We first compute the kinetic energy of a mode <j) n 
with amplitude A: Z(x,t) = A<j) n (x) cos ui n t. The speed 
of a mass element m c d.T is — A<f> n (x)uj n sin ui n t, the total 
kinetic energy is thus 



m c dx A 4> n (x)uj n sui oj n t 



2 



<f> n (x)dx I A u> n sin oj n t 



m c A 2 uj 2 sin' 1 u) n t 



(15) 
(16) 
(17) 



This is the kinetic energy of an harmonic oscillator of 
mass to c , resonant pulsation lo„ and amplitude A. No- 
tice that the amplitude of this oscillator is different 
from the deflection at the free end of the cantilever: 
A c = A4> n (x 1) = Az n (l) = Aa n ( n (l). We may 
also explicit the integral on </>^(x) in equation 16 using 
equation 14: 



1 



C, 2 {x)dx + toC^(1) A 2 w 2 sin 2 u> n t 



1 1 



—m c dx A 2 z' 1 l {x)uJ 2 l sin 2 uj n t 



-m b A c uj n sm u n t 



(18) 



(19) 



We easily identify here the sum of two terms: the kinetic 
energy of the mode n of the cantilever itself subject to a 
sinusoidal motion with an amplitude A c at its free end, 
and the kinetic energy of a point mass m;, subject to a 
sinusoidal motion with the same amplitude A c . The addi- 
tional Dirac term in (f> n thus takes into account the bead 
motion in the total energy of the equivalent harmonic os- 
cillator. The amplitude of the latter is not equal to the 
amplitude at the free end of the cantilever (which is also 
the case without the added mass since \4> n (x = 1)1 = 2 
for rh = 0). 

In the case where f ^ 0, no simple Dirac term can be 
added to the the eigenmode z n to ensure the orthogonal- 
ity of the basis. However, we can extend the energetic 
approach to include the kinetic energy due to the ro- 
tation of the bead, and ensure the normalization of the 
modes. Equation 14 is modified to the following criterium 
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FIG. 3: Thermal noise expected for a cantilever of static stiff- 
ness k = 0.1 N/m at 300 K for a few values of rh and f . The 
rms deflection of the first 5 modes is plotted as a function of 
the position along the cantilever x. The first mode is weakly 
impacted by the bead. When rh increases, the amplitude of 
the higher order modes vanish close to the free end, and the 
thermal energy is reported towards the antinodes, the ampli- 
tude of which raises. 



to compute the value of a n : 



C n {x)dx + rhQi{l) + fhf z a^) 



1 



(20) 



The result of this process is illustrated in figure 2 for a 
few values of rh and f. A similar approach has been used 
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by Oguamanam [19] to ensure the orthonormalization of 
the normal modes in a more general framework including 
the coupling of flexural and torsional modes. 



C. Thermal noise repartition between modes 

We compute the thermal noise of each resonant mode 
following [14, 20]: let us project the thermal noise driven 
deflection on the orthonormal basis <p n (x) : 



Z(x,t) = ^ Z n (t)<p n {x) 



(21) 



Under the hypothesis of uncoupled modes, we have for 
each degree of freedom 



h n (z 2 n (t)) = \k B T 



(22) 



where ks is Boltzmann's constant, T the temperature of 
the cantilever, and k n the stiffness of the mode defined 
by: 

k 

kn = -a* = m c u 2 n (23) 

The mean quadratic deflection measured in x should thus 
be 



(Z 2 (x,t))=^(Z„) 2 (t)|0„(x)| 5 



(24) 



ksT ^ \(j>n{x)\ 2 k B T ^ , 



n=l 



n=l 



(25) 



In figure 3, we plot the expected rms thermal noise at 
300 K along a cantilever of static stiffness k = 0.1 N/m 
for the five first eigenmodes, for a few values of in and 
f. Note that the normalization of the <j> n basis is a cru- 
cial step to apply the energy equipartition theorem in a 
quantitative manner to this analysis. Our approach also 
allows us to estimate the repartition of energy between 
the different modes. For example, if we perform the mea- 
surement at the free end of the cantilever (x = 1), the 
first mode account for r/i (x = l,m = 0,r = 0) = 97% 
of the total thermal fluctuations for fh = f = 0, and 
?7i(x = 1,771 = l,f = 0) = 99.8% for 771 = 1 and f = 0. 
In figure 4, we plot the contribution of each mode to 
the mean quadratic deflection measured at its extremity 
when 771 changes (for f = 0) : the larger the mass, the 
stronger is the contribution of the first mode. 



III. EXPERIMENTAL METHODOLOGY AND 
RESULTS 

A. Experiment description 

Manufacturer specifications of our two sample (A and 
B) are given in table I. Both present a "rectangular" 
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FIG. 4: Fraction r] n (x = 1, m,f = 0) at the free end of the 
cantilever of the mean quadratic deflection driven by thermal 
noise for the first 5 modes as a function of the normalized 
bead mass 771. The amplitude of the first mode (plotted as 
1 — 771) quickly overcomes the contribution of all the other 
modes, for which the inertia of the mass implies a node close 
to the free end of the cantilever. 




FIG. 5: Pictures from an optical microscope of beads glued 
at the apex of cantilevers A and B. 



geometry, close to the model used in our analytical ap- 
proach. However, the triangular end (see figure 5) de- 
parts from the model, and impedes a proper definition of 
their length L. We measure the deflection solely driven 
by thermal noise over the surface of these two differ- 
ent commercial cantilevers, first when they are still bare, 
then again after a glass bead has been glued at their free 
end. The radius r of the bead is reported in table I for 
each sample. 

As already mentioned, the uncertainty in the thickness 
of the lever is large due to the manufacturing process, re- 
sulting in a large uncertainty in the computation of its 
mass. The case is even worse for cantilever B, where the 
gold coating can change significantly the total mass due 
to the high density of gold. In addition, the quantity 
of glue cannot be measured precisely from the images of 
the cantilever. The geometric calculation of 777 can thus 
only give a rough estimation of the actual value. We 
estimate 777^4 rj 0.1 (range 0.04 — 0.3) for cantilever A, 
and 77773 w 1.2 (range 0.5 — 2.6) for cantilever B. The 
geometric estimation of f is less hampered by the uncer- 
tainty on the cantilever geometry, but still suffers from 
the uncontrolled repartition of the glue. We estimate 
r A = 0.02 ± 0.01 for cantilever A, and r B = 0.06 ± 0.02 
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Cantilever A 


Cantilever B 


Manufacturer 


Budget Sensors 


Nano- World 


Reference 


AIO, lever A 


Arrow-TL8Au 


Material 


Silicon 


Silicon 


Tip height 


(17 ±2) urn 


tipless 


Coating 


none 


Ti 5 ran + Au 30 nm 


Resonant frequency 


(15 ± 5) kHz 


6 kHz 
(3 - 14) kHz 


Force constant k 


0.2N/m 


0.03 N/m 




(0.04 - 0.7) N/m 


(0.004 - 0.54) N/m 


Length L 


(500 ± 10) Jim 


(500 ± 5) um 


Width W 


(30 ± 5) um 


(100 ± 5) um 


Thickness T 


(2.7 ± 1.0) um 


1.0 um 
(0.5 - 2.5) um 


Bead radius r 


(7.8 ± 0.2) jim 


(25.8 ±0.5) mri 



The complete set of results for cantil ever A is re ported 
in figure 6, where the rms amplitude y 7 (A^(x, y)) is rep- 
resented with a color coded scale. The first three vibra- 
tion modes can be clearly seen with their respective num- 
ber of nodes. A weak component in torsion can be seen 
for the third mode with the bead, certainly because it has 
not been glued perfectly on the axis. However, in the fol- 
lowing, we will neglect this effect and focus on the flexural 
modes along the x axis. Therefore, at each position x, 
the median yj (A^(x)) along the y axis of \J (A^(x, y)) is 
calculated. Due to the higher reflectivity of goald coated 
cantilever B, the background noise is lower and we accu- 
rately measure the thermal noise up to the fifth resonant 
mode for this last sample. 



TABLE I: Manufacturer specifications for cantilevers A and 
B. The last line corresponds to the glass bead that we glue at 
the free end of the cantilever, as illustrated in figure 5. 

for cantilever B. However, as demonstrated in [23], the 
tip of the bare cantilever A itself may have a none negli- 
gible effect. We can geometrically estimate »Mu P ~ 0.04 
and fAtip ~ 0.02 for a 17pm tall pyramidal tip. 

The measurement is performed with a home made in- 
terferometric deflection sensor [15], inspired by the origi- 
nal design of Schonenberger [21] with a quadrature phase 
detection technique [22] : the interference between the ref- 
erence laser beam reflecting on the chip of the cantilever 
and the sensing beam on the cantilever gives a direct mea- 
surement of the deflection with very high accuracy (see 
figure 1). This technique offers a very low intrinsic noise 
(down to 10~ 14 m/VHz). It is intrinsically calibrated as 
it measures directly the deflection against the wavelength 
of the laser beam, contrary to the standard optical lever 
technique that actually measures an angular deflection. 
At last, the focused beam size resolution is tuned to as 
small as 10 pm to ensure a good spatial resolution. 

We follow the same methodology as described in [15]. 
At every position x and y on a 5 pm x 5 pm grid, we mea- 
sure the deflection z(x, y, t) produced by the sole thermal 
excitation of the cantilever and we evaluate the power 
spectrum density (PSD) S z (x, y, f) on a 20s signal sam- 
pled at 2 MHz. For a quantitative characterization of the 
shape of the modes, the mean squared amplitude of each 
resonance (A^(a;,y)) is determined as a function of po- 
sitions x and y by integrating the PSD in a convenient 
frequency interval 2 A/ around each resonance frequency 

fn' 

rf n +Af 

(A 2 n (x,y))= / S 2 z (x,yJ)df. (26) 

This quantity is computed directly from the experimen- 
tal spectra, without any fitting process. We take care to 
subtract the background noise contribution of the inter- 
ferometer, and we also compensate for finite integration 
range in frequency [15]. 



B. Results and discussion 

1 . Resonant frequency ratios 

In a first attempt to estimate the added mass parame- 
ters m and f, which must be known for proper normaliza- 
tion of the normal modes (value of a n ), let us first focus 
on the ratio between the resonant frequency f n of the 
successive modes. Indeed, equation 23 translates into: 



a n {m,f) 
ai(m, f) 




The resonant frequencies of each mode are easily found 
by a simple harmonic oscillator fit of each resonant peak 
in the thermal noise spectrum. We report in table II 
the square root of the frequencies of the first 3 modes of 
cantilever A and the first 5 modes of cantilever B, nor- 
malized to the frequency of the mode 1. We can then 
compare those measurements to the output of the ana- 
lytical model (figure 11), and try to estimate the values 
of 77i and f for our samples. 

For cantilever A without the bead, the presence of the 
tip can be detected on the measured modes, and lead 
to the estimation m^tip = 0.057 ± 0.010 and f^tip = 
0.03 ± 0.02 (the standard deviations correspond to a 
0.2% max distance between the analytical model and 
measured frequency ratio on the 3 modes), close to the 
values expected from the geometrical analysis. Adding 
the bead change those values to vha = 0.336 ± 0.008 and 
¥a = 0.029 ± 0.005 (same criterium for error bars). 

For bare cantilever B, the ratios of frequencies between 
modes cannot be explained by our model : the triangular 
shape of the end of the cantilever alters the results cor- 
responding to a rectangular one. If however we look at 
the data for the loaded cantilever, we have a reasonable 
agreement for m B = 1.18 ± 0.09 and r B = 0.070 ± 0.006 
(the standard deviations correspond to a 2% max dis- 
tance between the analytical model and measured fre- 
quency ratio on the 5 modes), again in line with the ge- 
ometrical analysis. 
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FIG. 6: Maps of thermal noise of cantilever A. The rms amplitude of the first 3 flexural modes, with (right) and without (left) 
a bead glued at the free cantilever end, are plotted using a color coded amplitude map (color bar at the bottom of the figure). 
The resonant frequency and full scale amplitude are given below each map. A weak component in torsion can be clearly seen 
for the third mode with the bead, hinting at a slightly off axis gluing of the colloidal particle. 



2. Resonant frequency shifts - a naive attempt 



between successive modes of the loaded cantilever. 



Another way to determine the rh and f parameters is to 
analyse the resonant frequency shifts due to the loading, 
in an approach similar to the Cleveland method [24] . In 
a naive attempt, we suppose that the process of gluing 
a bead to the cantilever free-end should have a limited 
effect on its stiffness and proper mass, thus we can relate 
the frequencies with (f™) and without (f%°) the mass 
to the spatial eigenvalues with (a n (m, f)) and without 
(a„(0,0)) through equation 23: 



(rh, f) 



*n(0,0) 




(28) 



The resonant frequencies of each mode are easily found by 
a simple harmonic oscillator fit of each resonant peak of 
the thermal noise spectrum. We can then compare those 
measurements to the output of the analytical model, and 
try to estimate the values of rh and f for our samples. 

As shown in figure 10, the first mode is almost indepen- 
dent on the value of f, and should thus reliably be used 
to measure rh. This measurement is equivalent to the 
method of the added mass proposed by Cleveland [24]. 
We find with this protocol rh = 0.210 ± 0.005 for can- 
tilever A, and rh — 1.065 ± 0.010 for cantilever B. 

This estimation of rh is hardly compatible with the 
expectation of the higher order modes for cantilever A, 
which would rather be rh = 0.10 ±0.01 (mode 2) or 
rh = 0.055 ± 0.010 (mode 3), even considering the effect 
of r. As for cantilever B, the estimation could be com- 
patible with higher order modes, but for both cantilever 
the value of rh is clearly underestimated with respect to 
the previous measurement through the frequency ratios 



3. Resonant frequency shifts - refined analysis t 

Two naive hypotheses are responsible for the short- 
comings of the previous analysis of the frequency shifts 
upon loading: it first relies on the assumption that the 
effective length of the loaded cantilever is unchanged. It 
then assumes that the behavior of the cantilever without 
the load is that an ideal bare rectangular cantilever 

Let us first consider the effect of a possible modification 
of the effective length L of the cantilever upon gluing 
a bead close to its free end. Indeed, this process may 
rigidify the end portion of the cantilever, thus shortening 
its effective length by the rigid part. Alternatively, the 
inertia of a large mass not fixed exactly at the free end, 
by bringing the nodes of the higher order modes closer 
to the position of the bead than to the free end, can also 
lead to a shortening of the effective length of the loaded 
cantilever. In order to evidence the dependence of the 
eigenvalue a n in the cantilever's length L, equation 3 
can be expressed as: 



4 3//L 4 £< 

a„ = 



EI 



(29) 



where E is the cantilever's Young modulus, I its second 
moment of inertia, and \i its mass per unit length. As E, 
I and /j, do not depend on the cantilever's length nor on 
the gluing of a bead, equation 28 should in fact writes 



as 




(30) 
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Bare cantilever 


A: m A tip = 0.057, f A u P = 0.03 


B: m Bt ip = 0, fstip = 


Mode number n 


1 


2 


3 


1 


2 


3 


4 


5 


/■wo (Tl„\ 


14 956 


95 017 


267 530 


7195 


42 975 


117 220 


224 500 


363 370 


/ fwo / fwo 

Win Ih 


1.000 


2.521 


4.229 


1.000 


2.444 


4.036 


5.586 


7.107 


ainirhtip, nip)/Qi(m t ip, f t i p ) 
disagreement 


1.000 

0.0% 


2.521 
-0.0% 


4.236 
-0.1% 


1.000 
0.0% 


2.503 
-2.4% 


4.189 
-3.6% 


5.864 
-4.7% 


7.539 
-5.7% 


Loaded cantilever 


A: fh A = 0.35, f A = 0.03 


B: fh B = 1.18, r B = 0.06 


Mode number n 


1 


2 


3 


1 


2 


3 


4 


5 


Jn \ nL ) 


10 974 


83 050 


246 380 


3101 


29 319 


79 023 


152 400 


270 980 


V Jn III 


1.000 


2.751 


4.738 


1.000 


3.075 


5.048 


7.010 


9.348 


disagreement 


1.000 

0.0% 


2.762 
-0.4% 


4.762 
-0.5% 


1.000 
0.0% 


3.195 
-3.8% 


5.190 
-2.7% 


7.040 
-0.4% 


9.347 
0.0% 


Frequency shift upon loading 


cantilever A 


cantilever B 


Mode number n 


1 


2 


3 


1 


2 


3 


4 


5 


{L*/L*°)y/fS/Pr 


0.839 


0.916 


0.940 


0.643 


0.809 


0.805 


0.807 


0.846 


a n (m, f )/a n (m tip , f tip ) 

disagreement 


0.844 
-0.5% 


0.924 
-0.9% 


0.948 
-0.8% 


0.642 
0.2% 


0.820 
-1.2% 


0.796 
1.1% 


0.771 
4.7% 


0.796 
6.3% 



TABLE II: Frequency ratio between modes and frequency shift upon gluing of the bead. The square root of these ratios should 
be equal to the ratio of the corresponding spatial eigenvalues a n (fh, r) (with a correcting factor i w /i wo for the frequency shift, 
that we suppose equal to 0.98 here). Using estimated values of fh and r for each measurement set (bare and loaded cantilevers 
A and B), we get a good overall agreement for every mode, especially for sample A. The models reaches its limitation for the 
higher modes of cantilever B, whose triangular shaped end is not taken into account. 



where the superscript w (respectively wo ) designates a 
quantity with (respectively without) the load. Assuming 
oi%° is known, the parameters fh and f for the loaded 
cantilever can be deduced by comparing the a™ values 
obtained from equation 30 to the tabulated values given 
in appendix A [25]. Let us stress that if the resonant 
frequencies are measured with a very good accuracy, the 
effective length ratio is not known a priori. A small error 
in the effective length ratio leads only to the same rela- 
tive error on the value of a™. But since a n (fh,f) varies 
only slowly with fh (see appendix A), it can result in a 
much larger error on fh and thus on the normalization of 
the normalmodcs (parameter a n ). In the following, we 
estimate that the gluing of a bead reduces the effective 
length of cantilevers A and B by 10 pm, resulting in a 
2% shortening. We thus use this 2% correction in table 
II to compare the relative frequency shifts to the spacial 
eigenvalues ratios. 

When analysing the frequency ratios with and without 
a bead, it is also very important to use the right a™° 
value for the cantilever without the bead. We have seen 
with the naive approach that ignoring this initial load- 
ing for cantilever A leads to inconsistent values between 
modes, underestimating the true loading. However, if 
u n {mAt\ v ,r At - 1]? ) (with m At i p = 0.057 and f^tip = 0.03 
as determined in the first section) is used as the un- 
loaded reference value instead of cv„(0, 0), one obtains 
fh A = 0.37 ± 0.04 and f A = 0.031 ± 0.018 with the first 
three modes, in much better agreement with the previous 
estimation. 

The behavior of the tiplcss cantilever B also deviates 
from the one of a bare rectangular cantilever because 
of its triangular end. However, since the frequency ra- 



tios between the modes of the unloaded cantilever B do 
not yield any consistent set of fh and f values (see fig- 
ure 11), we choose to take mstip = ^stip = for the 
bare cantilever B. We thus compare the square root of 
the frequency ratio (corrected by the length ratio) with 
the a n (fh, f)/a n (0, 0) ratio displayed on figure 10 for the 
five first modes of cantilever B. As shown in figure 10, 
the first mode is almost independent on the value of f, 
and can thus be used to measure fh alone; we obtain 
ifiB = 1.17 ±0.04 for mode 1. The higher order modes 
are compatible with the estimation of itib = 1-17 and 
can be used to guess the value of fs- Using figure 10, 
we measure fs = 0.068 (mode 2), fs = 0.057 (mode 3), 
t b = 0.046 (mode 4), f B = 0.030 (mode 5). The disper- 
sion of results is quite large for fs-, and points to the lim- 
itations of the model with respect to the actual cantilever 
shape. A simultaneous least square minimization of the 
distance between the analytical model and measured fre- 
quency shifts on the 5 modes leads to : fhs = 1-19 ± 0.13 
and 7\b = 0.051 ± 0.006 (standard deviation correspond- 
ing to a 3.5 % max distance). 

As a summary, we have two ways to estimate the 
normalized added mass and equivalent gyration radius 
from the measurement of the resonant frequencies of the 
loaded and unloaded cantilever: the frequency ratio be- 
tween modes in one measurement, and the frequency 
shifts due to the addition of the bead. Provided the 
initial tip and effective length shortening are taken into 
account, both methods agree reasonably, though the dis- 
pertion on f is quite large. In the following we will retain 
the values : 

• Unloaded cantilever A: fh At i P = 0.057, f^tip = 
0.03. 
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FIG. 7: Amplitude of thermal noise for the first 3 flexural 
modes along the cantilever A with and without a sphere (in 
red and blue respectively). The markers represent the data 
while the lines exhibit the fits : independent fits of each mode 
in plain line, and simultaneous fit of all modes in dashed line. 
The agreement is excellent for the bare and the loaded can- 
tilever. 



• Loaded cantilever A: rfiA = 0.35, ?a — 0.03. 

• Unloaded cantilever B: metip = 0, retip = 0. 

• Loaded cantilever B: rag = 1.18, ¥b = 0.06. 

Eventually, using those estimations of fh and f , we can 
compute the values of the eigenvalues a n (rh,f) from the 
model, and compare them with the frequency shifts due 
to the bead. As one can read in table II, the agreement 
is quite good, with an overall agreement better than 1% 
for cantilever A, and 3% for cantilever B (except for the 
highest order modes of cantilever B, where the limitations 
of the model appear more severely). 



ft 



^ft 



TO 



100 




40 



20 





15 



10 




10 



o-i -r 



+ no bead 

indep. fit 

— — simult. fit 

o with bead 

indep. fit 

— — simult. fit 



























J 



















+ 










V 'r 















+ 








I ii 









t 








1 \ ■ 








i 

7 7 

/ 











100 200 300 400 500 
Position along the cantilever x (pi) 

FIG. 8: Amplitude of thermal noise for the first 5 flexural 
modes along the cantilever B with and without a sphere (in 
red and blue respectively). The markers represent the data 
while the lines exhibit the fits : independent fits of each mode 
in plain line, and simultaneous fit of all modes in dashed line. 
The agreement is very good for the bare and the loaded can- 
tilever. 



4- Spatial modes shapes 

Figures 7 and 8 display the rms amplitude y / (A^(x)) of 
the first 3 resonant modes of cantilever A and the first 5 
modes of cantilever B respectively, with (red circles) and 
without the bead (blue crosses) loading their end. It is 



worth mentioning that the maximum of the last mode is 
only a few pm high, demonstrating the high resolution of 
our instrument. The behavior of the experimental data 
is clearly in line with the model illustrated in figure 3: 
the effect of the bead is almost negligible on the first 
longitudinal mode, and the nodes of higher order modes 
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are shifted towards the free end of the cantilever. 

In a first step to fit the data, we use independent fits 
for all modes. The fitting function is equivalent to equa- 
tion 5. This generic function should be suitable for any 
cantilever clamped at its origin, whatever boundary con- 
ditions are applied at its other end. Note that x is now 
not normalized to the length L of the cantilever. We 
therefore use the following fitting function for (A^x)): 

ip {x,a,R, j-,x ) = 

f_x-x \ f_x-x 
a cos a — = — — cosh a 



-R 



X — Xq 



sinh a 



X — Xq 
L 



(31) 



with a, R, a/L and x~q the 4 fitting parameters. As can 
be seen in figures 7 and 8, the result of this procedure is 
excellent, the model closely matching the experimental 
data. 

The interesting output of the fit is the spatial eigen- 
value a n /L of each mode, however since L is not known 
precisely, only relatives values of a n can be compared to 
the theory. We choose to normalize the values of a n to 
the third mode: the presence of 2 nodes and the "high" 
amplitude constrain the fit to provide trustable values 
for CS3. The result of this procedure is displayed in tables 
III. The agreement with the theoretical ratios is good 
for higher order modes, within a few percent. Mode 1 
stands apart with a higher deviation: since no nodes are 
present, the fit is poorly constrained and the value of a.\ 
is not trustable. 



5. Stiffness determination 

If the generic independent fits are interesting to com- 
pare the eigenvalues and the shape of the cigenmodcs 
with the theory, it is not possible to compare the ampli- 
tude of the modes to the expectation from thermal noise 
excitation. Indeed, the normalization of the mode de- 
pends on the boundary conditions (kinetic energy of the 
bead), and cannot be guessed a priori. To go further, we 
therefore perform a simultaneous fit of all modes, impos- 
ing the values of a n [15]. The fitting function for (A^x)) 
arc now 



ipl(x,k,L,x ) 



a£(m,f) k 



L 



m,f ) (32) 



where 4> n {x, fn, f) are the normal mode defined and nor- 
malized in section II B (we write here explicitly the de- 
pendence in m and f to underline that those 4> n depend 
on the presence of the bead). The fitting parameters are 
k, L and afo, the values of fn and f are set to the estima- 
tion of section III B 3. 

We shall thus perform a single simultaneous fit of all 
modes with those 3 free parameters on each cantilever, 



bare or loaded with a bead. We tried several weighting 
of the modes to compensate the decreasing amplitude of 
higher order modes driven by thermal noise: the function 
6j to minimize during the fit is defined as 



(k,L,x ) =^a& / dx\(A 2 n (x)) - <p 2 n (x, k, L, x )\ 



(33) 

where j is a weighting parameter: j — corresponds to 
natural weighting (mode 1 dominant), and increasing j 
weights more and more the higher order modes ("flat" 
weights for j = 8). We estimate the best fitting parame- 
ters k, L and x~o for j = to 12, to test the robustness of 
the simultaneous fitting procedure. The dashed lines in 
figures 7 and 8 represent the result of this fitting process. 
Though not as perfect as independent fits, the results are 
in good agreement for all modes of each cantilever, with 
and without the bead. 

For the bare cantilever A, the best fit values are k = 
(0.321 ± 0.008) N/m and L = (497 ± 1) urn. The uncer- 
tainties correspond to the standard deviation in the full 
range of weighting parameter j, their low values demon- 
strating the robustness of the fit. For the loaded can- 
tilever A, the best fit values are k = (0.339 ± 0.010) N/m 
and L = (485 ± 7) pm. The stiffness of the cantilever 
experiences a small increase (6%) after the gluing of the 
bead. This is not surprising since the glue increases the 
rigidity of the end of the cantilever, shortening its effec- 
tive length by about 10 pm according to the fit. This 2% 
decrease in length translates into a 6% rise of the stiffness 
(k scales as 1/L 3 ), in agreement with our estimation. 

For the bare cantilever B, the best fit values are 
k = (0.151 ± 0.002) N/m and L = (463 ± 4) pm. Again, 
the dispersion of estimated parameters is very low, hint- 
ing at the robustness of the model and fitting proce- 
dure. The value of L is quite small with respect to 
the manufacturer specifications, however this effective 
length takes into account the triangular end of the can- 
tilever. For the loaded cantilever B the best fit values are 
k = (0.169 ± 0.010) N/m and L = (448 ± 5) pm. Again, 
this 12% increase of the static spring constant goes in 
the expected direction, and is coherent with a reduction 
of the effective length of the cantilever by 15 pm (which 
should translate into a 10% increase in k). 



IV. CONCLUSION 

Our work demonstrates that even a glued bead as large 
as 10% of the length of a soft cantilever (k ~ 0.1 N/m) 
modifies only slightly its first flcxural mode and its static 
stiffness. A simultaneous fit of the thermal modes lead 
to a small decrease in the effective length and a small in- 
crease in the effective static stiffness of a cantilever upon 
gluing the bead. However, since these variations respect 
the 1/L 3 scaling of k, the stiffness at the geometrical free 
end of the cantilever appears to be the same with and 
without a bead loading. Thus, the classic method to de- 
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Bare cantilever 


A: friAtip = 0.056, fUtip = 0.03 


B: 




0, fstip = 


= 


Mode number n 


1 


2 


i 
i 


2 


4 


5 


Measurement: cv n /c*3 


0.278 


0.599 


n 077 
U.z / { 


0.602 


1.389 


1.760 


Theory: a n (m t - lp , f t ip)/a3(m t ip, ftip) 


0.236 


0.595 




0.598 


1.400 


1.800 


disagreement 


-18% 


-0.7% 


-1070 


-0.7% 


0.8% 


2.2% 


Loaded cantilever 


A: fh A = 0.35, f A = 0.03 


B: 


rh B = 1.18, f B = 0.06 


Mode number n 


1 


2 


1 


2 


4 


5 


Measurement: a n /c\3 


0.251 


0.580 


0.191 


0.561 


1.389 


1.865 


Theory: a n (m,f)/a3(rh,f) 


0.210 


0.580 


0.193 


0.616 


1.356 


1.801 


disagreement 


-19% 


0.0% 


-1.1% 


8.8% 


-2.7% 


-3.5% 



TABLE III: Eigenvalues normalized to the third mode for both cantilevers, bare and loaded with the bead: we compare the 
output of generic fits with equation 31 and the values expected from the model. The agreement is quite good for every modes 
but the first one, where the fit is not sufficiently constrained to converge to accurate estimations. 



tcrmine the spring constant by measuring the thermal 
spectrum of the fisrt fiexural mode at the free end of the 
cantilever can still be used even if the microlever is func- 
tionalized. In fact, we have shown that the mode 1 based 
thermal calibration gets even closer to the static stiffness 
when the size of the bead increases, since the first mode 
gathers a higher fraction of the thermal energy at the free 
end. We have also shown that the mass model consider- 
ing the load as a modification of the boundary conditions 
at the free end of a beam is a good approximation and 
fit well all our results. On the way, we have introduced 
a proper normalization method of the resonant modes, 
an unavoidable step to compute the thermal noise am- 
plitude of each mode. Once the stiffness of the colloidal 
probe determined through the classic thermal noise cal- 
ibration, one will make sure to take into account other 
corrections due to the tip geometry, as demonstrated in 
reference [26]. 

One of the main difficulties in testing the models is to 
estimate correctly the mass and gyration radius of the 
glued bead. We have seen that it can be important to 
consider the initial loading due to the AFM tip and the 
effective length shortening in order to get consistent re- 
sults: an error of 30% on fh could have been made by 
only considering the frequency shift due to the bead, if 
one refers to a tiplcss cantilever. The ratio between reso- 
nant frequencies of a cantilever offers an interesting way 
to estimate the load properties (without any prior knowl- 
edge on the cantilever, like its unloaded properties), by 
comparison to tabulated values of a n (fh,f). We provide 
a set of such values in appendix A, for the first 5 modes, 
< fh < 2 and < f < 0.1 

As a final remark, let us study how our findings ap- 
ply to the common AFM detection scheme. Indeed, our 
differential interferometer allows us to measure the ac- 
tual deflection in any point of the cantilever, whereas 
most AFM use an optical angular deflection measure- 
ment. The latter technique needs to be calibrated to 
infer from the 4 quadrants photodiode output the true 
deflection of the probe. This step is usually done by a 
rigid contact between the probe and a hard surface, the 
calibrated displacement of the sample providing a bench- 
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Normalized tip mass fh 



FIG. 9: Multiplicative correction factor to apply to the ther- 
mal noise measurement of the first resonant mode with an an- 
gular deflection measurement technique to compute the can- 
tilever stiffness: the higher the added mass, the smaller the 
correction. 



mark to measure the sensitivity (in nm/V). This sensi- 
tivity is thus valid for a static deformation only, and a 
mode dependent correction factor must be applied to es- 
timate the actual deflection of a resonant mode [1]. This 
multiplicative factor is 0.817 when one wants to use the 
thermal noise measurement of mode 1 to calibrate the 
spring constant of a classic cantilever. In the current 
framework, we can easily compute how this correction 
factor depends on the normalized bead mass fh by com- 
paring the slope of the first eigenmode to that of a static 
deflection with the same deflection at its free end. We 
plot the result in figure 9 (computed with r = 0): the 
correction rapidly vanishes as fh increases. Indeed, we 
have seen that the effect of the added mass is to decrease 
the eigenvalue a\ (see figure 10), thus equation 2 tends to 

= when fh increases. This last equation is that de- 
scribing the static deflection, thus the first normal mode 
tends to the static deflection when fh increases. The ef- 
fect of the added mass is thus twice in the good direction 
for the angular measurement: both the sensitivity of the 
sensor tends to be more accurate (figure 9), and the first 
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mode gathers most of the thermal noise at the cantilever 
free end (figure 4). More generally, an accurate coeffi- 
cient can be extracted from our analysis for any rh, and 
applied to the thermal noise calibration of AFM colloidal 
probes in any commercial devices [27]. 
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Appendix A: Spatial eigenvalues a„(rh,r) 

In this appendix, we plot in figure 10 the spatial eigen- 
values a n (rh,r) numerically computed for modes 1 to 5 



as a function of rh and f , normalized to their value at 
rh = f = 0. The curve for a% is very close to the pre- 
diction given in the Cleveland method [24], linking the 
frequency shift to the added mass. This figure can be 
used to estimate rh and f from the value of the frequency 
shift of various modes, if the initial situation corresponds 
to a rectangular tipless cantilever and the potential effec- 
tive length decrease after gluing the particle is known. 

In figure 11, we plot the same computed eigenvalues 
a n (m, f), but normalized to the value of the first mode 
for the same added mass ct\ {rh, f) . When one has no prior 
knowledge of the unloaded resonant frequencies, the ratio 
of the resonant frequencies between modes can be used 
to estimate rh and f. 

Finally, tables IV to VIII report numerical values of 
a n (rh,f), for the first 5 modes, < rh < 2 and < f < 
0.1. 
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FIG. 10: Spatial eigenvalues a TO (m,f) numerically computed 
for modes 1 to 5 as a function of rh and f , normalized to their 
value at rh = 0. The values of these ratios for cantilevers A 
and B, estimated by the frequency shift due to the addition 
of the bead and supposing an effective length decrease of 2%, 
are plotted as a labelled horizontal line. Mode 1 is almost 
independent in the value of f, and can thus be used to estimate 
rh. This curve is very close to the prediction given in the 
Cleveland method [24]. 
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FIG. 11: Spatial eigenvalues a„(m,f) numerically computed 
for modes 1 to 5 as a function of m and f, normalized to the 
value of mode 1. The values of these ratios for cantilevers A 
and B, estimated by the frequency ratios between modes, are 
plotted as a labelled horizontal line for unloaded (index wo ) 
and loaded (index w ) cantilevers. The model is not suited 
to the description of bare cantilever B (B wo ), but allows the 
estimation of rh and f for the 3 other situations. 
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Mode 1: ai(m, f) 
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TABLE IV: Qi(m, f): table of eigenvalues of mode 1 for < m < 2 and < r < 0.1 . 



Mode 2: 012 (m, r) 
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9980 


3 


9669 


3 


9290 


3 


8846 


3 


8343 


3 


7788 


3 


7192 







90 


4 


0411 


4 


0373 


4 


0260 


4 


0071 


3 


9806 


3 


9465 


3 


9049 


3 


8564 


3 


8016 


3 


7416 


3 


6775 


jfl 


1 


00 


4 


0311 


4 


0271 


4 


0148 


3 


9942 


3 


9653 


3 


9282 


3 


8831 


3 


8305 


3 


7714 


3 


7070 


3 


6386 




1 


10 


4 


0228 


4 


0184 


4 


0051 


3 


9829 


3 


9517 


3 


9116 


3 


8629 


3 


8064 


3 


7431 


3 


6745 


3 


6022 




1 


20 


4 


0157 


4 


0109 


3 


9967 


3 


9728 


3 


9393 


3 


8962 


3 


8440 


3 


7836 


3 


7162 


3 


6437 


3 


5677 




1 


30 


4 


0096 


4 


0045 


3 


9892 


3 


9637 


3 


9278 


3 


8818 


3 


8262 


3 


7619 


3 


6907 


3 


6144 


3 


5350 




1 


40 


4 


0042 


3 


9988 


3 


9826 


3 


9554 


3 


9172 


3 


8683 


3 


8092 


3 


7412 


3 


6662 


3 


5863 


3 


5038 




1 


50 


3 


9995 


3 


9938 


3 


9766 


3 


9477 


3 


9073 


3 


8554 


3 


7929 


3 


7213 


3 


6427 


3 


5594 


3 


4740 




1 


60 


3 


9954 


3 


9893 


3 


9711 


3 


9406 


3 


8978 


3 


8430 


3 


7772 


3 


7021 


3 


6200 


3 


5335 


3 


4454 




1 


70 


3 


9916 


3 


9853 


3 


9661 


3 


9340 


3 


8889 


3 


8312 


3 


7620 


3 


6834 


3 


5980 


3 


5086 


3 


4179 




1 


80 


3 


9883 


3 


9816 


3 


9615 


3 


9277 


3 


8803 


3 


8197 


3 


7473 


3 


6653 


3 


5767 


3 


4845 


3 


3915 




1 


90 


3 


9853 


3 


9783 


3 


9572 


3 


9217 


3 


8720 


3 


8086 


3 


7330 


3 


6477 


3 


5560 


3 


4611 


3 


3660 




2 


00 


3 


9826 


3 


9752 


3 


9531 


3 


9161 


3 


8640 


3 


7978 


3 


7190 


3 


6306 


3 


5359 


3 


4385 


3 


3414 



TABLE V: a2{m, f): table of eigenvalues of mode 2 for < m < 2 and < f < 0.1 . 



Mode 3: as(m, f) 
f 
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5756 
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8913 


6 


7015 


6 


4737 


6 


2386 


6 


0212 


5 


8331 


5 


6762 


5 


5477 




1 


30 


7 


1199 


7 


0940 


7 


0126 


6 


8689 


6 


6653 


6 


4251 


6 


1830 


5 


9637 


5 


7772 


5 


6237 


5 


4994 




1 


40 
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8470 
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7262 
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4562 
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50 
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1134 
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0838 
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9905 
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8255 


6 


5954 


6 


3339 


6 


0814 


5 
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5 


6794 


5 


5333 
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60 


7 


1108 


7 


0793 


6 


9800 


6 


8043 


6 


5616 


6 


2910 


6 


0349 


5 


8153 


5 


6364 
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4941 
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3820 
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70 
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1084 


7 


0751 
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9697 


6 


7834 


6 


5286 


6 


2498 


5 


9910 


5 


7726 


5 


5968 


5 


4583 


5 
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80 
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1063 


7 


0712 
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9597 
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7627 


6 


4964 


6 


2103 
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9496 
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7328 


5 


5602 
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4255 
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3208 
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90 
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1044 
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0675 
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7423 
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4649 
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5 


9104 
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3953 
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2941 
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7 
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4341 
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TABLE VI: a 3 (m, f): table of eigenvalues of mode 3 for < fa < 2 and < f < 0.1 . 



Mode 4: 04 (fa, f) 
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10.9955 


10.9955 


10.9955 


10.9955 


10.9955 


10.9955 


10.9955 







10 


10 


5218 


10 


5127 


10.4849 


10.4372 


10.3676 


10.2745 


10.1581 


10.0210 


9.8698 


9 


7129 


9.5591 







20 


10 


4016 


10 


3862 


10.3387 


10.2555 


10.1325 


9.9701 


9.7771 


9.5714 


9.3720 


9 


1921 


9.0377 







30 


10 


3480 


10 


3269 


10.2607 


10.1427 


9.9681 


9.7450 


9.5000 


9.2643 


9.0586 


8 


8894 


8.7547 







40 


10 


3178 


10 


2910 


10.2061 


10.0527 


9.8278 


9.5550 


9.2796 


9.0376 


8.8423 


8 


6913 


8.5764 







50 


10 


2984 


10 


2660 


10.1620 


9.9726 


9.7003 


9.3899 


9.1010 


8.8654 


8.6860 


8 


5530 


8.4547 







60 


10 


2850 


10 


2469 


10.1235 


9.8976 


9.5823 


9.2460 


8.9550 


8.7317 


8.5690 


8 


4519 


8.3669 







70 


10 


2751 


10 


2313 


10.0882 


9.8257 


9.4728 


9.1205 


8.8348 


8.6260 


8.4788 


8 


3751 


8.3008 







80 


10 


2675 


10 


2181 


10.0547 


9.7563 


9.3715 


9.0111 


8.7349 


8.5409 


8.4077 


8 


3152 


8.2494 







90 


10 


2615 


10 


2064 


10.0225 


9.6889 


9.2779 


8.9156 


8.6511 


8.4714 


8.3503 


8 


2672 


8.2083 


rii 


1 


00 


10 


2566 


10 


1958 


9.9911 


9.6236 


9.1919 


8.8321 


8.5803 


8.4138 


8.3033 


8 


2280 


8.1749 




1 


10 


10 


2526 


10 


1860 


9.9602 


9.5603 


9.1128 


8.7587 


8.5199 


8.3654 


8.2641 


8 


1954 


8.1471 




1 


20 


10 


2492 


10 


1768 


9.9297 


9.4992 


9.0403 


8.6941 


8.4680 


8.3243 


8.2309 


8 


1679 


8.1237 




1 


30 


10 


2463 


10 


1681 


9.8994 


9.4403 


8.9738 


8.6369 


8.4229 


8.2890 


8.2026 


8 


1444 


8.1037 




1 


40 


10 


2438 


10 


1597 


9.8694 


9.3837 


8.9128 


8.5861 


8.3836 


8.2585 


8.1781 


8 


1242 


8.0864 




1 


50 


10 


2417 


10 


1516 


9.8396 


9.3293 


8.8568 


8.5408 


8.3490 


8.2317 


8.1568 


8 


1065 


8.0713 




1 


60 


10 


2398 


10 


1438 


9.8099 


9.2773 


8.8054 


8.5002 


8.3184 


8.2082 


8.1380 


8 


0910 


8.0580 




1 


70 


10 


2381 


10 


1361 


9.7804 


9.2275 


8.7580 


8.4636 


8.2912 


8.1874 


8.1214 


8 


0772 


8.0463 




1 


80 


10 


2366 


10 


1286 


9.7511 


9.1799 


8.7145 


8.4307 


8.2668 


8.1687 


8.1066 


8 


0650 


8.0358 




1 


90 


10 


2352 


10 


1212 


9.7220 


9.1345 


8.6743 


8.4008 


8.2449 


8.1521 


8.0933 


8 


0540 


8.0264 




2 


00 


10 


2340 


10 


1139 


9.6931 


9.0913 


8.6371 


8.3737 


8.2251 


8.1370 


8.0813 


8 


0440 


8.0179 



TABLE VII: ai(fa,f): table of eigenvalues of mode 4 for < fn < 2 and < f < 0.1 . 
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Mode 5: as(m, f) 
f 

O00 O01 O02 003 004 O05 006 O07 008 009 OTTO 






00 


14 


1372 


14 


1372 


14 


1372 


14 


1372 


14 


1372 


14 


1372 


14 


1372 


14 


1372 


14 


1372 


14 


1372 


11 


1372 





10 


13 


6142 


13 


5953 


13 


5364 


13 


4314 


13 


2743 


13 


0670 


12 


8284 


12 


5890 


12 


3745 


12 


1966 


12 


0553 





20 


13 


5067 


13 


4742 


13 


3700 


13 


1781 


12 


8973 


12 


5727 


12 


2730 


12 


0358 


11 


8619 


11 


7373 


11 


6475 





30 


13 


4615 


13 


4160 


13 


2659 


12 


9859 


12 


6053 


12 


2337 


11 


9480 


11 


7515 


11 


6196 


11 


5296 


11 


4665 





40 


13 


4367 


13 


3782 


13 


1803 


12 


8138 


12 


3660 


11 


9938 


11 


7428 


11 


5833 


11 


4805 


11 


4116 


11 


3637 





50 


13 


4210 


13 


3494 


13 


1019 


12 


6552 


12 


1717 


11 


8212 


11 


6053 


11 


4741 


11 


3911 


11 


3359 


11 


2975 





60 


13 


4102 


13 


3253 


13 


0267 


12 


5097 


12 


0149 


11 


6941 


11 


5082 


11 


3981 


11 


3291 


11 


2833 


11 


2514 





70 


13 


4023 


13 


3040 


12 


9532 


12 


3778 


11 


8882 


11 


5980 


11 


4367 


11 


3425 


11 


2837 


11 


2447 


11 


2175 





80 


13 


3963 


13 


2843 


12 


8808 


12 


2594 


11 


7854 


11 


5236 


11 


3821 


11 


3002 


11 


2491 


11 


2151 


11 


1914 





90 


13 


3916 


13 


2658 


12 


8096 


12 


1539 


11 


7012 


11 


4646 


11 


3392 


11 


2669 


11 


2218 


11 


1918 


11 


1709 


1 


00 


13 


3878 


13 


2479 


12 


7398 


12 


0603 


11 


6314 


11 


4169 


11 


3047 


11 


2401 


11 


1998 


11 


1730 


11 


1542 


1 


10 


13 


3846 


13 


2305 


12 


6715 


11 


9774 


11 


5730 


11 


3777 


11 


2764 


11 


2181 


11 


1817 


11 


1574 


11 


1404 


1 


20 


13 


3820 


13 


2134 


12 


6051 


11 


9040 


11 


5237 


11 


3449 


11 


2527 


11 


1997 


11 


1665 


11 


1444 


11 


1288 


1 


30 


13 


3797 


13 


1966 


12 


5408 


11 


8390 


11 


4816 


11 


3172 


11 


2327 


11 


1841 


11 


1536 


11 


1333 


11 


1190 


1 


40 


13 


3778 


13 


1798 


12 


4788 


11 


7811 


11 


4454 


11 


2934 


11 


2155 


11 


1707 


11 


1425 


11 


1237 


11 


1104 


1 


50 


13 


3761 


13 


1632 


12 


4193 


11 


7296 


11 


4139 


11 


2728 


11 


2007 


11 


1590 


11 


1329 


11 


1154 


11 


1030 


1 


60 


13 


3746 


13 


1465 


12 


3622 


11 


6836 


11 


3863 


11 


2549 


11 


1876 


11 


1489 


11 


1244 


11 


1081 


11 


0965 


1 


70 


13 


3733 


13 


1299 


12 


3078 


11 


6424 


11 


3620 


11 


2391 


11 


1762 


11 


1398 


11 


1170 


11 


1016 


11 


0907 


1 


80 


13 


3721 


13 


1132 


12 


2558 


11 


6052 


11 


3404 


11 


2251 


11 


1660 


11 


1318 


11 


1103 


11 


0958 


11 


0856 


1 


90 


13 


3711 


13 


0966 


12 


2065 


11 


5718 


11 


3212 


11 


2125 


11 


1569 


11 


1247 


11 


1043 


11 


0906 


11 


0810 


2 


00 


13 


3701 


13 


0798 


12 


1596 


11 


5414 


11 


3039 


11 


2013 


11 


1487 


11 


1182 


11 


0989 


11 


0860 


11 


0768 



TABLE VIII: a 5 (m,f): table of eigenvalues of mode 5 for < m < 2 and < f < 0.1 . 



